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ABSTRACT 

We present a comprehensive study of the observational constraints on spatially flat cosmological models 
containing a mixture of matter and quintessence — a time varying, spatially inhomogeneous component 
of the energy density of the universe with negative pressure. Our study also includes the limiting case of 
a cosmological constant. We classify the observational constraints by red shift: low red shift constraints 
include the Hubble parameter, baryon fraction, cluster abundance, the age of the universe, bulk velocity 
and the shape of the mass power spectrum; intermediate red shift constraints are due to probes of 
the red shift-luminosity distance based on type la supernovae, gravitational lensing, the Lyman-alpha 
forest, and the evolution of large scale structure; high red shift constraints are based on measurements 
of the cosmic microwave background temperature anisotropy. Mindful of systematic errors, we adopt 
a conservative approach in applying these observational constraints. We determine that the range of 
quintessence models in which the ratio of the matter density to the critical density is 0.2 <J D, m <; 0.5 
and the effective, density-averaged equation-of-state is — 1 < w <; —0.2, are consistent with the most 
reliable, current low red shift and microwave background observations at the 2a level. Factoring in the 
constraint due to the recent measurements of type la supernovae, the range for the equation-of-state 
is reduced to — 1 < w < —0.4, where this range represents models consistent with each observational 
constraint at the 2a level or better (concordance analysis). A combined maximum likelihood analysis 
suggests a smaller range, — 1 < w <J —0.6. We find that the best-fit and best-motivated quintessence 
models lie near fi m w 0.33, h « 0.65, and spectral index n s — 1, with an effective equation-of-state 
w w —0.65 for "tracker" quintessence and w = — 1 for "creeper" quintessence. 



1. INTRODUCTION 

The most widely studied cosmological models at the 
present time are variants of the Cold Dark Matter (CDM) 
paradigm within which adiabatic perturbations in a dom- 
inant CDM species grow due to gravitational instability 
from quantum fluctuations imprinted during an inflation- 
ary era. The bulk of the evidence today strongly favors 
models within which Q m < 1 and any hot com ponent is 
not significant, Ohom *C Ocnivr- Several authors (Ostrikcrl 
& Steinhardt 1995|, fCrauss & Turner 1995|, [Turner & Whit 



problem and prediction of a nearly scale-invariant spec- 
trum of energy density fluctuations. Finding these ar- 
guments compelling , we will adopt the ansatz that the 
universe is spatially flat. 

In our previous work ( Ostriker fc Steinhardt 1995| ) we 



1997[ ) have found that the best and simplest fit concordant 



with current observations is provided by 

w [0.30 ±0.10] + [0.04 ±0.01] 



(1) 



with h = 0.65 ± 0.15. One is thus led to either an open 
universe or one in which the remaining energy density re- 
quired to produce a geometrically flat universe is some 
additional energy component (E) with 



i. 



(2) 



An important advantage of the flat model is that it is con- 
sistent with standard inflationary cosmology, and its asso- 
ciated resolution of the cosmological horizon and flatness 



identified the range of models consistent with then-current 
observations, restricting attention to the case where £Ie is 
vacuum energy or cosmological constant (A). The cos- 
mological constant, a static, homogeneous energy compo- 
nent with positive energy density but negative pressure, 
was introduced initially by Einstein in a flawed attempt 
to model our universe as a static spacetime with positive 
spatial curvature. In a spatially flat universe, however, the 
negative pressure results in a repulsive gravitational effect 
which accelerates the expansion of the universe. This ear- 
lier analysis, which preceded by several years the recent 
evidence based on luminosity-red shift surveys of type la 
supernovae, already found that A is favored over standard 
cold dark matter or open models. Over the intervening 
years, the supernovae measurements as well as other ob- 
servations have reinforced these conclusions. 

In this paper, we update our earlier concordance anal- 
ysis and expand it to include the possibility that the ad- 
ditional energy component consists of quintessence, a dy- 
namical, spatia lly inhomogeneous fo rm of energy with neg- 
ative pressure (Caldwell et al. 1998). A common example 
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Lindc 1999 and references 
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is the energy of a slowly evolving scalar field with positive 
potential en ergy, similar to the inflaton field in inflationary 
cos mology flWeiss 19871 iRatra k Peebles 198j jWcttcricfj 



1995 1, iFricman ct al. 19951, Coble ct al. 19971, Ferrcira k 

Joyce [l997|, |Caldwcll ct al. 1998|, ^latev et al. 1998|) . Unlike al. 1998| that there is a one-to-one correspondence between 



size that scalar field models of quintessence are not only 
the simplest, well- motivated choice from a particle physics 
standpoint, but, also, they can mimic fluids wi th arbitrary 
equation -of-state. (It was shown explicitly in Caldwell ct 



a cosmological constant, the dynamical field can support 
long-wavelength fluctuations which leave an imprint on 
the cosmic microwave background (CMB) and large-scale 
structure. In particular, the long wavelength fluctuations 
change the relation between the amplitude of the CMB 
anisotropy and the gravitational potential fluctuations so 
that the COBE normalization of the mass power spectrum 
depends on the pressure of the quintessence component. A 
further distinction is that w, the ratio of pressure (p) to 
energy density (p), is — 1 < w < for quintessence whereas 
w is precisely —1 for a cosmological constant. Hence, the 
expansion history of the universe for a A model versus a 
quint essence model (for the same 0„, today, say) is diffcr- 



ent. In general, the acceleration, the age and the volume 
of the universe arc less for quintessence models than for A 
mode js" (assuming all other cosmic parameters are fixed). 



A prime motivation for considering quintessence models 
is to address the "coincidence problem," the issue of ex- 
plaining the initial conditions necessary to yield the near- 
coincidence of the densities of matter and the quintessence 
component today. For the case of a cosmological con- 
stant, the only possible option is to finely tune the ratio 
of vacuum density to matter-radiation density to 1 part 
in ~ 10 120 at the close of inflation in order to have the 
correct ratio today. Symmetry arguments from particle 
physics are sometimes invoked t o explain wh y the cosmo- 
logical constant should be zero (Banks 1996) but there is 
no known explanation for a positive, observable vacuum 
density. For quintessence, because it couples directly to 
other forms of energy, one can envisage the possibility of 
interactions which may cause the quintessence component 
to naturally adjust itself to be comparable to the mat- 
ter densi ty today. In fact, recen t investigations (Zlatev et 



al. 19^, Steinhardt et al. 1999] ) have introduced the no 



tion of "tracker field" m odels of quintessence which have 
attra ctor-like solutions (Peebles k Ratra 1988, Ratra k 
PeebL p~1988 ) which produce the current quintessence en- 
ergy density without the fine-tuning of initial conditions. 
A related development h as been "creeper field" models 
(Huey k Steinhardt 199E) which are nearly as insensitive 
to initial conditions but indistinguishable from A today. 

Fundamental physics provides some further motivation 
for light scalar fields. Particle physics theories with dy- 
namical symmetry breaking or non-perturbative effects 
have been found which generate potentials with ultra- 
light ma sses which support negative pressure (Affleck ct 
al. 19 §5l |Hill k Ross 1988a|, |Hill k Ross 1998b|, | Binctruy 



et al. 1 19961 , [Barrciro ct al. 1998[ [Binctruy 1999| ). These 
suggestive results lend appeal to a particle physics basis 
for quintessence, as a logical alternative to an ad hoc invo- 
cation of a cosmological constant. We do not aim to base 
our investigation of the properties of quintessence cosmolo- 
gies on a specific particle physics model, however, as such 
models are still in a developmental stage. An intriguing 
thought is that progress in the cosmological observations 
and experiments discussed here will soon decide the is- 
sue, possibly pointing to new fundamental physics inac- 
cessible in the accelerator laboratory. We would empha- 



a general time-dependent equation-of-state and an equiva- 
lent scalar field potential.) Ge neralizati ons to tensor fields 
or mo re general stress tens or (Hu 1999) or topological de- 
fects ( Spcrgcl &; Pen 1996| ) have also been considered; for 



the purposes of this study based on current observations 
they can be well described by scalar fields, in addition. 

Extending the realm of cosmological models to in- 
clude quintessence opens up a new degree of freedom, 
the equation-of-state of the missing energy component. 
The added degree of freedom necessarily makes the pro- 
cedure of selecting viable models more complicated. Pre- 
vious studies ha ve considered some spe c ific combinations 
of ob s ervations ( Silvcira k Waga 1997 Turner fc White 



1997| , Parnavich et al 1998|, |Perlmutter et al 1998iror 



some specifi c models (Ferreira fc Joyce 1997, Ferrcira k 



Joyce 1998). Here we systematically examine the most 



complete range yet of measures of astrophysical phenom- 
ena at low, intermediate, and high red shift and make a 
complete search in parameter space to objectively identify 
the viable models. We show that the observations are con- 
sistent with A and quintessence models for a substantial 
range of parameters. An impressive feature is how a num- 
ber of observations, not only the measurements of type la 
supernovae, favor a missing energy component with sub- 
stantially negative pressure. 

The pace of cosmological observations is proceeding so 
rapidly that any quantitative conclusions may soon be- 
come dated. Nevertheless, we think that an assessment at 
the present time is worthwhile for at least three reasons. 
First, our study shows a concordance among a growing 
numb er of observations. Compare d to the previous anal- 
ysis (Ostriker k Steinhardt 1995) new constraints have 
been added and old ones have been revised, and, yet, no- 
tably, the key conclusions — a new energy component and 
an accelerating expansion rate — have been significantly 
strengthened. Second, the study isolates those observa- 
tions which are playing the lead roles in shaping the cur- 
rent conclusions and identifies observations or combina- 
tions thereof which will be most decisive in the near fu- 
ture. Third, this comprehensive analysis enables one to 
identify specific best-fit models which can be explored in 
much greater detail to search for more subtle implications 
and tests. 

The organization of the paper is as follows. In section ^ 
we discuss the parameterization of the cosmological con- 
stant and quintessence cosmological models. In section [| 
we present the observational constraints, classified by low, 
intermediate, and high red shift. We evaluate the con- 
straints, presenting the results in section [|. We conclude 
in section ^| with an identification of the overall best-fit 
models and a discussion of future observations. 

2. PARAMETERIZATION OF QUINTESSENCE 
COSMOLOGICAL MODELS 

The quintessence (QCDM) cosmological scenario is a 
spatially-flat FRW space-time dominated by radiation at 
early times, and cold dark matter (CDM) and quintessence 
(Q) at late times. For simplicity, we will consider models in 
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which the quintessence component consists of a scalar field 
slowly rolling down its effective potential with a constant 
equation-of- state. The detailed equations-of-motion are 
discussed in Caldwell et al. 1998. This class of models is a 



good approximation for most of the range of quintessence 
candidates. The models we consider can then be fully 
characterized by the following five parameters: 

wq: A constant quintessence equation-of-state, in the 
range w G [— 1, 0] in the present epoch. In most cases, the 
quintessence equation-of-state changes slowly with time, 
but the observational predictions are well approximated 
by treating w as a constant, equal to 



dafi,Q(a)w(a)/ / daQ,Q(a) 



(3) 



We will comment on exceptions later. Tracker field 
quintessence models have a lower bound on the range of 
w, while creeper quintessence leads to w very close to —1, 
which is effectively indistinguishable from a true cosmo- 
logical constant. 

fi m : The matter density parameter, defined as the ratio 
of the matter energy density, including CDM and baryons, 
to the critical energy density p c = 3H 2 /8irG, where H is 
the Hubble constant. We assume that any contribution to 
the energy density due to a hot component, such as neu- 
trinos, is small, insofar as free-streaming has a negligible 
effect on the clustering of CDM. Unless otherwise specified, 
we have imposed J7 m + Qq = 1, where Qq is the corre- 
sponding density parameter for quintessence. Hence, the 
matter density parameter lies in the range tt m £ [fib, 1]- 

fib'. The baryon density parameter, defined as the ra- 
tio of the baryonic energy density to the critical energy 
density. 

h: The Hubble parameter, related to the Hubble con- 
stant by H = 100 h km/s/Mpc. 

n s : The index of the power spectrum of primordial den- 
sity fluctuations in the matter and radiation. This param- 
eter also controls the contribution of tensor perturbations, 
for which we consider two cases. In the first case, we im- 
pose the inflationary relation between the amplitude of the 
primordial density and gravitational wave (tensor) pertur- 
bations for n s < 1, revised for the case of quintessence 
( Caldwell fc Steinhardt 1998| ). For n s > 1, we assume the 
tensor contribution is negligible. In the second case, we 
assume the tensor contribution is negligible for all values 
of the spectral index. (We only illustrate the first case; the 
second case yields an indistinguishable result in regard to 
our concordance analysis.) 

The most revealing way to depict the concordance of 
constraints on quintessence models is to project the five- 
dimensional parameter space into the £l m -w plane. In dis- 
playing this plane, we assume the universe is spatially flat. 
Since the flatness condition requires = 1 — £7 m , the pa- 
rameters w and fl m completely specify the quintessence 
portion of the cosmology. ACDM corresponds to the line 
w = — 1, and SCDM corresponds to the line f2 m = 1 in 
this plane. 

3. OBSERVATIONAL CONSTRAINTS 

We take a conservative approach in applying the cosmo- 
logical constraints. Observational cosmology is currently 
in a period of rapid growth so that the current constraints 



must be considered as work in progress, rather than final. 
Certainly, measurements of many astrophysical phenom- 
ena are becoming more refined, with greater precision as 
statistical errors are reduced, and with greater accuracy 
as systematic errors are better understood. If this situa- 
tion described all observations, we could confidently apply 
the results to the full limits of the published errors. Be- 
cause the observational constraints typically restrict com- 
binations of independent model parameters, then by com- 
bining several experimental results, we could find tighter 
parameter bounds than if the observations were applied 
individually. For example, combining two constraints re- 
stricts a two-parameter system to an ellipse in parameter 
space, whereas the individual constraints applied succes- 
sively allow a rectangle which contains that ellipse. The 
former combination would allow the determination of a 
model which was "best" in the maximum likelihood sense. 
However, not all measurements have well controlled sys- 
tematic errors or even high precision statistical errors. Nor 
do they have true gaussian errors. Nor are the errors un- 
corrected. A prime example is the current set of type 
la supernovae magnitude-red shift data, which must be 
interpreted with caution. For these reasons, we do not ad- 
vocate cosmological parameter extraction with the current 
set of astrophysical data. Until the coming generation of 
precision experiments are on-line, we believe it can be mis- 
leading to combine the full set of astrophysical data as if 
the errors were statistical and gaussian. 

For this reason, we employ an additional procedure, 
which we call "concordance," to evaluate the observational 
constraints. We identify models as passing the observa- 
tional constraints if they lie within 2er of each individual 
constraint. (We do not consider the joint probabilities 
spanning two or more observations.) We allow a generous 
range for systematic errors. Not only does this procedure 
provide a reliable picture of current constraints, but there 
is the added advantage that it renders transparent which 
observations are most important in delimiting the range 
of currently allowed parameters and which future observa- 
tions will be most influential. 

One might characterize the difference between the two 
approaches as follows: two identical and observationally 
independent constraints on a single parameter do not 
change, at all, the allowed range of that parameter in the 
concordance analysis, but do reduce the range in a max- 
imum likelihood approach. Clearly, the concordance ap- 
proach is too conservative if the errors are known to be 
gaussian. However, in the current case, systematics domi- 
nate and the correlations between errors in different obser- 
vations is unclear. In these circumstances, as discussed in 
the Appendix, the maximum likelihood approach can pro- 
duce seriously misleading results. One should be especially 
watchful and examine closely situations where the concor- 
dance and maximum likelihood analyses strongly disagree. 
Maximum likelihood has the undesireable feature that it 
will seek a compromise among conflicting data. Hence, 
we advocate the more conservative concordance approach, 
supplemented with comparison to a full maximum likeli- 
hood estimator for the constraints assuming gaussian er- 
rors. 

We classify the observations by red shift. At low red 
shift, the constraints are due to: the Hubble con- 

stant; age of the universe; baryon density; x-ray cluster 
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abun dance; shape of the mass power spectrum; and bulk (Buries fc Tytler 1997a, Buries fe Tytler 1997b, Buries & 



flow. A_t intermediate red shift. 



1, the constraints are 



due to: type la supernovae; evolution of the x-ray cluster 
abundance; gravitational lcnsing; Lyman-a forest. At high 
red shift, z»l, the constraints are due to: the fluctuation 
amplitude and spectral tilt based on the large angle CMB 
temperature anisotropy (COBE); small angle CMB tem- 
perature anisotropy. The high red shift CMB constraint 
due to COBE, along with the set of low red shift results, 
serve as the strongest, most reliable constraints. These 



ytler 1998)^ If this value reflects the primordial abun- 
dan ce, then big bang nucleos ynthesis (BBN; for a review 
see [Schramm fc Turner 199§| and references therein) with 



three light neutrinos gives ilt,h 2 = 0.019 ± 0.002, where 
the la error bars allow for possible systematic uncer- 
tainty. While other observations suggest that the abun- 
dance varies on galactic length scales where it is expected 
to be uniform, suggesting that heretofore unknown pro- 
cesse s may be processing the deuterium ( Jenkins ct alj 



consti aints will be shown to dominate the boundary of 1999), we will adopt the hypothesis that the cosmological 



the allowed parameter range. The intermediate red shift 
constraints, due to SNe and the evolution of x-ray clus- 
ters, are rapidly reaching the point where they impact the 
range of cosmological models. The small angle CMB mea- 
surements, as well, are soon to yield prime cosmological 
information. We will consider each of these constraints in 
turn. 

The observational constraints used to restrict the 
quintessence parameter space are listed in the following 
subsection. The low red shift, and COBE-based high red 



abundance has been ascertained by the measurements of 
Buries & Tytler. 

BF: Observations of the gas in clusters have been used 
to estimate the baryon fraction (compared to the total 
mass) t o be /"„„., = (0.06 ± 0.003)fe~ 3/2 flEvrard 1997 ; 
also see [White et al. 1993j [Fukugita et al. 1997|) . The 



shift c onstraints compose the core set of concordance tests 
of our cosmological models. The remaining intermedi- 
ate and high red shift constraints are less certain at the 
present, although they offer the promise of powerful dis- 
crimination between models in the near future. 

3.1. Low Redshift 

H: The Hubble constant has been measured through 
numerous techniques over the years. Although there has 
been a marked increase in the precision of extragalactic 
distance measurements, the accurate determination of H 



stellar fraction is estimated to be less than 20% of the 
gas fraction, so that f stellar = 0.2/i 3 / 2 / gQS . Next, sim- 
ulations suggest that the baryon fraction in clus ters is 
less than t he cosmological value by about 10% ( Lubir 



et al. 1996 ) representing a depletion in the abundance of 
baryons in clusters by a factor of 0.9 ±0.1. Hence, the 
cosmological baryon fraction (f2i/f2 m ) is estimated to be 
f baryon = (0.067 ± 0.008)/i~ 3/2 + 0.013 at the la level. Us- 
ing the observed baryon density from BBN, we obtain a 
constraint on f2 m : 



0.019/1" 



0.067/i- 3 / 2 + 0.013 



(l±0.32) 



(4) 



at the 2a level. For h = 0.65, this corresponds to a value 
of n m = 0.32 ±0.1. 

The baryon fraction has also be en estimated on smaller 
scales such as galaxy systems (see McGaugh 1998 and ref- 
erences therein). However, the relationship between the lo- 
cal baryon fraction on those scales and the global baryon 
fraction is uncertain and currently beyond the power of 
numerical study. 

a&: The abundance of x-ray clusters at z = provides 
a model dependent normalization of the mass power spec- 
trum at the canonical 8/i _1 Mpc scale. The interpretation 
fe Kochanek 1997| , |Kundic ct al. 1997| , |Schcchtcr ct alj of x-ray cluster data for the ca se of quintessence models 
1997| ). Using surface brightness fluctuations to calibrate has been carried out in detail in |Wang fe Steinhardt 1998 



has been slow. The Ho Key Project ( Freedmanet al. 1998| ), 
which aims to measure the Hubble constant to an accu- 
racy of 10%, currently finds H — 73 ± 6(stat) ± 8(sys) 
km/s/Mpc; the method of type la supernovae gives H = 
63.1 ±3 4(internal) ± 2.9(external) km/s/Mpc ( |Hamuy ct| 
al. 1996); a recent measurement based on the Sunyaev- 
Zcldovich ef fect in four nearby clusters gives H = 54 ± 14 
km/s/Mpc (Myers et al. 1997); typical values obtained 
from gravitational lens systems arc H ~ 50 — 70 km/s/Mpc 
with up to ~ 30% errors ([Falco ct al. 19971- I KcctorJ 



the bright cluster galaxy Hubble diagram, a far-field value 
H = 89 ± 10 km/s/Mpc, ou t to ~ 11, 000 km/s, has been 
obtained (Lauer et al. 1998). Clearly, convergence has not 
yet been reached, although some methods are more prone 
to systematic uncertainties. Based on these diverse mea- 
sures, our conservative estimate for the Hubble parameter 
is H = 65±15 km/s/Mpc at effectively the 2a level. While 
knowledge of H would certainly be a decisive constraint, 
our current uncertainty will not prevent us from narrowing 
the field of viable cosmological models. 

t : Recent progress in the dating of globular clusters 
and the calibration of the cosmic distance ladder has re- 
laxed the lower bound on the age of th e universe . We 



in which case the constraint is expressed as 

a 8 W m = (0.5 -0.19) ±0.1 
where the error bars are 2a, with 



7 : 

e 



: 0.21 - 0.22k; - 
: (n s -l) + (h- 



0.33O-, 
-0.65). 



0.250 



(5) 



(0) 



adop t tp > 9.5 Gyr as a 95% lower limit (Chaboyer et al 



1998 , galaris fc Weiss 1998 ), although we note that some 
arguments ( |Paczynski 199S ) suggest a 10% higher limit is 
more appropriate. 

BBN: Recent observations of the deuterium abundance 
by Buries and Tytler yield D/H = 3.4 ± 0.3(stat) x 10~ 5 



This fitting formula is valid for the range of parameters 
considered in this paper. 

Perhaps the two most important constraints on the mass 
power spectrum at this time are the COBE limit on large 
scale power and the cluster abundance constraint which 
fixes the power on 8/i -1 Mpc scales. Together, they fix 
the spectral index and leave little room to adjust the power 
spectrum to satisfy other tests. 

Shape: If light traces mass with a constant bias factor 
on large scales, then the deprojected APM galaxy cluster 
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data ( Peacock 1997b ) can be used to constrain the shape 
of the underlying mass power spectrum. The bias factor 
is defined as b 2 = Papm /Pmass, the ratio of the APM to 
mass power spectra on a given scale; we assume the bias 
represents a constant, quadratic amplification of the clus- 
tering power of rare objects over the density field on large 
scales ( Kaiser 1986 ). Hence, in keeping with our spirit 
of conservativism, we restrict our attention to wavenum- 
bers which are well within the linear reg ime, using only the 
seven lowest frequency bins (as given in Peacock 1997a| ) for 
scales above 8 Mpc/h. The shape constraint consists of the 
requirement that the mass power spectrum fit the seven 
APM data points with b > 1 and a reduced \ 2 5= 2.0, corre- 
sponding to a confidence level of 95%. In effect, our shape 
test depends also on the power spectrum amplitude. The 
lower bound on the bias is due to the assumption that the 
brigh t, luminous APM objects are preferentially formed 



in hig hly ovcrdcnsc regions (jPavis et al 1985|, | Bardee 



1986 



tively. Extrapolating to smaller scales, within the class 
of ACDM models, the SFI constraint can be recast as 
<j$f = 0.82 ±0.12. Due to the discrepancy with the cluster 
abundance constraint, we are hesitant to apply this recent 
result until further analysis reinforces its conclusions. 

The bulk flow on the largest scales provides another 
method to do cosmology with the velocity field. In Figure|| 
we compare the predictions of a set of QCDM models with 
observation. The large sample variance on the Maxwellian 
distributed velocity field means that consistency requires 
the observations lie below the upper 95%CL bulk velocity. 
(The lower 95%CL bulk velocity is very small, so we may 
effectively treat this constraint as an upper bound.) A 
measurement near or above the swath of predicted curves 
could serve as a strong indicator of the cosmology. For 
comparison, we a lso show the measure d bu lk velocities of 
|Dckcl ct al. 1999 [, piovanelli ct al. 1998 , and Lauer & Post 



Kaiser 1986, Cen & Ostriker 1998, Blanton et al 



1998) (Although the bias may be very large at the time of 
formation, simple arguments indicate that by the present 
time, b may have decreased to no lower than unity. See 



Fry 1986, Tegmark fc Peebles 1998| ). The consequences of 
b < 1 will be discussed. While we give no upper bound on 
6, almost all the best fitting values for concordance mod- 
els fall within the rough (model-dependent) upper bound 
estimate of b < 1.5 based on higher-order statistics of the 
APM data set and current theoretical ab initio modeling 
( Gaztanaga fc Fricman 1994 , |Ccn fc Ostriker 1998 , Blan- 



ton et 



al. 



& Wjllick 
1997" 



1 



man 1994. At present, the diversity of measurements, as 
displayed in Figure |2[ dilutes the strength of the constraint 
resulting from the comparison of the velocity dipole with 
the CMB dipole. 

3.2. Intermediate Redshift 

SNe: Type la supernovae are not standard candles, 
but empirical calibration of the light curve - luminos- 
ity relationship suggests that the objects can be used 
as distance indicators. There has been much progress 
in these observations recently, and there promises to be 
more. Hence, a definitive constraint based on these re- 
sults would be premature. However, we examine the re- 
cent results of th e High- Z Supernova Sea rch Team (HZS: 
Riess ct al. 1998 , Garnavich et al. 1998) and the Sup er- 
nova Cosmology Project (SCP: Pcrlmutter et al. 1998| ) to 
constrain the luminosity distance - red shift relationship in 
quintessence cosmological models. We have adopted the 
following data analysis procedure: we use the supernova 
data for the shape of the luminosity - red shift relationship 
only, allowing the calibration, and therefore the Hubble 
constant, to float; we excise all SNe at z < 0.02 to avoid 
possible systematics due to local voids and overdensities; 
for SNe at z > 0.02, we assume a further uncertainty, 
added in quadrature, corresponding to a peculiar velocity 
of 300km/s in order to devalue nearby SNe relative to the 
more distant ones (for the SCP data, a velocity of 300km/s 
has already been included). There is substantial scatter in 
the supernovae data, as seen in Figure ||. The scatter is 
1993|, IWillick et al. 1997| , |Willick fc Strauss- so wide that no model we have tested passes a \ 2 test 

with the full SCP data set; using a reduc ed set, Fit C, 
argued b y the SCP as being more reliable Pcrlmutter et] 



1998). Our computations show that the popular 
"shape parameter" T = VL m h is not an accurate description 
of the goodness of fit to the APM data, given the variety 
of models that we consider here, since the amplitude of the 
power spectrum is characterized by other combinations of 
parameters, including w. As illustrated in Figure [ll five 
sample models with T ranging from 0.20 to 0.52 all pass 
our shape test based on a % 2 analysis. 

Velocity Field: The large-scale velocity field has long 
been used as a means to probe the background and fluc- 
tuation matter density field. One method is to compare 
peculiar velocity data obtained from distance indicators, 
such as Tully-Fisher, and from red shift surveys, in or- 
der to estimate il m , modulo an assumption about biasing. 
Through this method, the quantity (3 = f(fl m )/b is ob- 
tained, where / = dhi5 m /dhia ~ fl 1 ^, 6 and b is a linear 



bias parameter. A variety of recent results ([Davis, Nusser 



a Costa et al. 1998) find (3 ~ 0.5 - 0.6. For a bias 



not too different from unity, these results suggest a low 
matter density. In the lack of a more complete understand- 
ing of bias, however, this method cannot be transformed 
into a rigorous constraint on J7 m . As well, it remains to 
be understood why similar approaches which compare the 
density fields obtained from distance indicators and red 
shift survey data, along the lines of the POTENT method 
( Bcrtschinger fc Dekel 198S), genera lly obtain higher val- 
ues, e.g. f3 ~ 0.9 ( |3igad ct al. 1998 ). Another method is 
to compare observations with the predicted velocity field 
within the context of a particul ar cosmological model . Re- 
sults based on the Mark III (Zaroubi ct al. 1997| ) and 
SFI QZehavi 1998| , Freudling ct al 1999| ) catalogs yield 
the constraint f z P(k) = (4.8 ± 1.5) x 10 3 (Mpc/h) 3 and 
(4.4 ± 1.7) x 10 3 (Mpc/h) 3 at k = O.lh/Mpc, respec- 



al. 1998, a finite range of models do pass the x z test and 
the range is comparable to the range obtained by the \ 2 
test using the HZS data set. To gauge the current situa- 
tion, we will report both \ 2 tests and maximum likelihood 
tests; to be conservative, we use the largest boundary (the 
X 2 test based on HZS data using MLCS analysis) for our 
concordance constraint. 

Cluster Evolution: The abundance of rich clusters — 
objects presumed to have formed from high density peaks 
drawn from the exponential tail of an initially Gaussian 
perturbation distribution - can be used to constrain the 
amplitude of the mass power spectrum at intermediate red 
shift. The current observations have been converted into a 
number density of clusters above a certain mass threshold 
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Mi. 5, defined to be the mass within the comoving radius 
Rcom = l-5/i -1 Mpc. For the models of interest, the abun- 
dance evolves approximately as a power law for < z < 1: 



n(> Ml 5 ,z) oc 10 A(Ml - : 



(7) 



see |Wang fe Steinhardt 1998|) . The bigger A{M X , 5 ) is, the 



weaker the evolution is, implying lo w Q m and w. Since 



the measurements (summarized in Bahcall et al. 1997 ) 
are still in the preliminary stage, we adopt A{M\,5 = 
8 x 10 14 ft, _1 Mpc) = -l.Ttxi as a conservative, 2a limit. 
Similar tests hav e been applied in the context of ACDM 
and open models (Bahcall et al. 1997, Carlberget al. 1997), 
for which the results arc model dependent. 

Lensing Counts: The statistics of multiply imaged 
quasars, lensed by intervening galaxies or clusters, can be 
used to determine the luminosity distance - red shift rela- 
tionship, and thereby constrain quintessence cosmological 
models. There exists a long literature of est imates of the 
lensing constraint on A m odels (e.g. from Turner et al 



1984[To Falco et al. 1998). In one approach, the cumula- 
tive lensing probability for a sample of quasars is used to 
estimate the expected number of lenses and distribution 
of angular separations. Using the Hubble Space Te lescope 
Snapshot Survey quasar sample ( Maoz et al. 1993), which 
found four lenses in 502 sources, Maoz- Kix ( {Maoz fc Rix 
19930 arrived at the limit n A £ 0.7 at the 95% CL. In 
a series of stu dies, similar con s traints have bee n obtained 
using o ptical (Kochanck 1995| , Kochanck 1996) and radio 
lenses (|Falco et al. 1998|). Waga and c ollaborators ( Tor-| 
res & Waga 1996 , |Waga fc Miceli 199S ) have generalized 
these results, finding that the constraint weakens for larger 
values of the b ackground equation-of- state, w > — 1 (as 
noted earlier by Ratra & Quillen 1992). In our evaluation 
of the constraint based on the HST-SSS data set, we find 
that the 95% confidence level region is approximately de- 
scribed by J7q <J 0.75 + (1 + w) 2 , until the inequality is 

Torres fc Wag^ 



satu rated at w = —1/2, c onsistent with 



1996, Waga fc Miceli 1999. In principle, t 



lis test is a sen- 



tive best-fit models discussed in Table I. The mass power 
spectra in the upper panel all satisfy COBE normaliza- 
tion at large scales and the cluster abundance constraint 
on as on 8 hr 1 Mpc scales. Since the models have already 
passed the constraints on 8 hr 1 Mpc scales and higher, 
one might hope that tests of the linear power spectrum at 
smaller scales might further distinguish the models. How- 
ever, smaller scales correspond to the non-linear growth 
regime where effects like scale-dependent bias make it dif- 
ficult to compare observations to the linear power spec- 
trum. Measurements of the Lyman-a forest are promising 
because they probe the power spectrum on smaller scales 
at a red shift before non-linearities develop and, hence, en- 
able direct comparison to the linear power spectrum. We 
draw the reader's attention to the fact that, in converting 
to z = 2.5, we have made a model-dependent rescaling of 
the abscissa so that the units are those of velocity, which 
then allows direct comparison to the data. In the upper 
(z = 0) panel, the models differ substantially on large 
scale but appear to converge on small scales. Projecting 
back to z = 2.5 and rescaling, one might hope that the 
models are distinct due to the differing growth functions. 
However, instead, examples with the same spectral tilt 
and n m Ho/ 'H(z) (as is demonstrated inadvertently by our 
representative models) nearly overlap everywhere, making 
discrimination very difficult. When fl m HQ/H(z) is fixed, 
as shown in the figure, then, we note, the Lyman-a mea- 
surements can be used to determine the tilt, n s . 

The current data appears to favor n s = 1. There is some 
hope that improved limits can discriminate among mod- 
els with different tilt, but determining other parameters, 
and especially discriminating between cosmological con- 
stant and quintessence on this basis, appears hopeless be- 
cause the predictions of various models converge, as shown 
in the lower panel. Of course, this conclusion applies not 
only to Lyman-a forest measurements, but any approach 
that measures the mass power spectrum at moderate red 
shift. 

Further applications of the Ly-a fore st, such as the 



sitive probe of the cosmology: however, it is susceptible to abun dance of damped Ly-a absorbers ( pardner et al 



a nu mber of systematic errors (for a discussion, see Mal- 
hotra let al. 1997 . Cheng fc Krauss 199S ). Uncertainties 



1997) and correlations in the lines of sight at red shifts 



in the j luminosity function for source and lens, lens evo- 
lution, lensing cross section, and dust extinction for opti- 
cal lenses, threaten to render the constraints compatible 
with or even favor a low density universe over fl m = 1. 
Taking the above into consideration, none of the present 
constraints on quintessence due to the statistics of multi- 
ply imaged quasars are prohibitive: models in concordance 
with the low-z constraints are compatible with the lensing 
constraints. 

Ly-a (and other mass power spectrum measurements 
at moderate red shift, 1 < z < 10): The Ly-a forest 
has b een used as the basis of a number of cosmological 



2-4 (McDonald fc Miralda-Escude 1998, Hui et al 



1999) have been developed as tests of geometry and cx- 



probes. Most recently, the effect of the local mass den- 



pansion history, although no substantial constraints have 
as yet been obtained. 

There are a number of observational probes which are 
sensitive to the cosmology, but which have not yet ma- 
tured into critical tests. We list some of these tests which 
may prove to be powerful constraints in the near future. 
Measures of the abundance of objects, similar to the clus- 
ter evolution constraint, can be used to gauge the growth 
of stru cture. Observation s of galaxies formed as early as 
z ;> 3 ( Stcidcl et al. 1998| ) have been interpreted, on the 
basis of a Press-Schechter formalism (Press fc Schechtei 



sity in the i ntergalactic mediu m on the Ly-a optical depth 
(Hui 1999, Croft et al. 1998a) has been used to estimate 
the m ass power spectrum at a r ed shift of z = 2.5 (Croft et 
al. 19)8b, [Weinberg et al. 199£ ). This is a good pedagogi- 
cal ex ample to study what can and cannot be learned from 



1974 ), to suggest that among a family of CDM cosmolo- 
gies, flat, low-density models best satisfy the constraint 



studies at moderate red shift 1 < z < 10. In Figure 4 we 
show the linear mass power spectrum today (upper panel) 
and at red shift z = 2.5 (lower panel) for the representa- 



(|Mo fc Fukugita 1996|) . 

Finally, it has been proposed to use the statistics of grav- 
itational lens arcs produced by intermediate red shift c lus- 
ters as a me a ns to distinguish cosmo logical models ( Wu 
fc Mao 1996 , Bartelmann et al. 1998 ). These are poten- 



tially extremely powerful tests, insofar as the rare, high 
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density fluctuations reflect the underlying cosmology In 
particular, the lens arcs statistics are exponentially sen- 
sitive to the growth function, which differs for A versus 
quintessence models and, hence, has the potential of dis- 
tinguishing the two scenarios. However, results based on 
numerical simulations are not yet capable of resolving the 
core of clusters accurately enough to provide reliable limits 
( Wambsganss et al. 1998) from the theoretical side. 



3.3. High Redshift 

One of the most powerful cosmological probes is the 
CMB anisotropy, an imprint of the recombination epoch 
on the celestial sphere. The large angle temperature 
anisotropy pattern recorded by COBE can be used to place 
two constraints on cosmological models. 

COBE norm: The observed amplitude of the CMB 
power spectrum is used to constrain the amplitude of the 
u nderlying density per turbations. We adopt the method 



explain the observed anisotropy spectrum. At present, 
however, there is ample CMB data which can be used to 
constrain cosmological models. 

We take a conservative approach in applying the small 
angular scale CMB data as a model constraint. Our inten- 
tion is to simply determine which quintessence models are 
consistent with the ensemble of CMB experiments, rather 
than to determine the most likely or best fitting model; ex- 
amining Figure ^|, the error bars are clearly so large that 
a 'best fit' has little significance. We have restricted our 
attention to that subset of published CMB experimental 
data which satisfies the following objective criteria: multi- 
frequency; positive cross-correlation with another experi- 
ment; careful treatment of foregrounds. (While wc advo- 
cate these criteria, our conclusions are not strongly sen- 
sitive to this selection of data.) Hence, we use the band- 



of Bu]|m fc White 19"9"7 to normalize the power spectrum to 



COBl^C^^eTise^^Scunca version oi CMBEAST (|ScT^ 

fc Zalctarriaga 19^o[ ) to compute the CMB anisotropy spec- 
tra, this normalization is carried out automatically.) We 
have verified that this method, originally developed for A 
and open CDM models, can be applied to the quintessenc e 
cosmological models considered in this work (Dave 1998). 
Of course, there is uncertainty associated with the COBE 
"normalization" : the 2er uncertainty in r ms quantities is 
appro ximately 20% (see footnote #4 in jBunn fc White 
1997), which conservatively allows for statistical errors, as 
well as the systematic uncertainty associated with the dif- 
ferences in the galactic and ecliptic frame COBE map pix- 
elizations, and potential contamin ation by high-latitude 
foregrounds(e.(7. Gorski et al. 1998] ). 

n s : COBE has been found to be consistent with a 



power estimates fro m COBE ( Bennett et al. 1996| ), Python 
([Piatt et al. 1997|) MSAM (Ichcng ct al. 19*971), QMAP 
(|Pcvlin ct al. 1998^ |Herbig et al. 1998|, |01ivcira-Costa~el 
al. 1998;), S askatoon iNetterfield et al. 1997|) CAT ([Scott 



= 1.2 ± 0.3 spectral index (Gorski et al. 1996, Hin 



shaw |et al. 199TS J, but this assumes the only large an 
guiar scale anisotropy is generated via the Sachs- Wolfe 
effect on the last scattering surface. This neglects the 
baryon-photon acoustic oscillations, which produce a rise 
in the spectrum, slightly tilting the spectrum observed by 
COBE. In general, the spectral index determined by fitting 
the large angular scale CMB anisotropy of a quintessence 
model, which is also modified by a late-time integrated 
effect, to the shape of the spectrum tends to overesti- 
mate the spect ral tilt. For example , analysis of a class of 
CDM models ( [Hancock et al 1998j ) (ACDM and SCDM, 
a subset of the models considered here) finds a spectral 
tilt n s = 1.1 ± 0.1. We conservatively restrict the spec- 
tral index of the primordial adiabatic density perturba- 
tion spectrum, with P(k) oc k n " , to lie in the interval 
n s s [0.8,1.2]. Note that inflation generically predicts 
n s ~ 1, with n s slightly less than unity preferred by infla- 
ton potentials which naturally exit inflation. 

Small Angle CMB: Dramatic advances in cosmology 
are expected in the near future, when the MAP and Planck 
satellites return high resolution maps of the CMB tem- 
perature and polarization anisotropy. When the measure- 
ments are analyzed, we can expect that the best deter- 
mined cosmological quantities will be the high multipole 
Ci moments, such that any proposed theory must first 

5 www.cita.utoronto.ca/ ~knox/radical/bpdata.html 
6 www. sns.ias.edu/~max/cmb/experiments. html 



ct al. 1996D, and RLNG5M (|Leitch et al. 1998[) experiments 



as the basis of the small angle CMB cosmological con- 
straint. Figure ^| shows the bandpower averages at the ef- 
fective multipole number, l e , with several QCDM models 
for comparison. We apply a simple \ 2 test with the pre - 
dicted bandpower averages, STi c (defined in Bond 1995). 
(A compilation of bandpower averages and window func- 
tions is available from Knox 5 and Tegmark 6 .) Since the 
reported bandpower errors are typically not Gaussian dis- 
tributed, treating the STi c as a Gaussian random variable 
can introduce a bias in the estimation of the quality of 
agreement. In this case, we also consider a y 2 test in the 
quantity ln(£TfJ, following pond et al. 1998j . 

Figure |^ illustrates features relevant to limits on the 
spectral index, n s . The COBE li mit on the spectral i ndex 
is reported to be n s — 1.2 ± 0.3 (Gorski ct al. 1996, Hin-| 



shaw et al. 1996) based on comparison between data and 



standard CDM models. In the Figure, though, it is appar- 
ent that there is negligible difference in the low angular 
scale predictions despite a range of An = 0.2. The dif- 
ference in the large angular scale integrated Sachs- Wolfe 
contribution compensates for the difference in spectral in- 
dex. Hence, as will be addressed in a later paper, the 
COBE limit on n s is somewhat expanded when QCDM is 
included. 

4. CONCORDANCE RESULTS 

We have evaluated the cosmological constraints for the 
set of quintessence models occupying the five dimensional 
parameter space: w, fi m , fif,, h, n s . The results are best 
represented by projecting the viable models onto the 
Q m — h and Q m — w planes. 

The concordance region due to the suite of low red shift 
constraints, including the COBE normalization and tilt 
n s , are displayed in Figures |], [?]. Each point in the shaded 
region represents at least one model in the remaining three 
dimensional parameter space which satisfies the observa- 
tional constraints. 

In Figure ^, the boundaries in the Q, m direction are de- 
termined by the combined BBN and BF constraints as a 
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function of h, while h is only restricted by our conserva- 
tive allowed range and the age constraint. The age does 
not impact the f2 m — h concordance region, since for the 
allowed values of f2 m and h, there is always a model with 
a sufficiently negative value of w to satisfy the age con- 
straint. Relaxing either the BBN or BF constraint would 
raise the upper limit on the matter density parameter to 
allow larger values of f2 m . This requires a simultaneous re- 
duction in the spectral index, n s , in order to satisfy both 
the COBE normalization and cluster abundance. 

In Figure [?], the upper and lower bounds on fl m are 
again determined by the combination of BBN, BF, and 
h. The lower bound on Q m due to the combination of 
the BBN and BF constraints can be relaxed if we allow 
a more conservative ran ge for the baryon density, such as 



0.006 < fl b h 2 < 0.022 ( |Levshakov 1998| , |Qlive, Stcigmarj 



& Wa lker 1999| ). However, the constraints due to ag and 
the shape of the mass power spectrum take up the slack, 
and the lower boundary of the concordance region is rel- 
atively unaffected. The lower bound on 17 m near w = — 1 
is determined in part by the shape test; the mass power 
spectrum in a model with low fl m and strongly negative 
w is a poor fit to the shape of the APM data, based on a 
X 2 -test. This constraint on models near w = —1 is relaxed 
if we allow anti-bias (b < 1), although b < 1 is strongly 
disfavored on a theoretical basis. At the other end, for 
w ;> —0.6, the lower bound on £l m is determined by the 
combination of the upper bound on the spectral index, and 
the x-ray cluster abundance constraint on as- If we further 
restrict the bias to b < 1.5, a small group of models at the 
upper right corner with w J> —0.2 and fl m ;> 0.4 will fail 
the shape test. 

We see that models occupying the fraction of the param- 
eter space in the range — 1 < w <; —0.2 and 0.2 <; fl m <; 
0.5 are in concordance with the basic suite of observations, 
suggesting a low density universe. It is important to note 
that the set of viable models spans a wide range in w; the 
concordance region is not clustered around w — —1, or 
A, but allow such diverse behavior as w ~ —1/3. How- 
ever, the case w = 0, which can result from the scaling 
expo nential potential ( Ratra fe Peebles 1988) , Peebles & 
Ratra 1988), is clearly in contradiction with observation: 
the fi m required by the x-ray cluster abundance constraint 
is incompatible with the matter density parameter allowed 
by th e BF and BBN constraints. Hence, the models with 



C explored in Fcrreira & Joyce 1997, Ferreira & Joyce Miceli 199S , as well as the mor e sophisticated analyses car 



1998 are not viable. 



the constraint dominating a particular portion of the con- 
tour. However, we show that by lifting the shape test, the 
constraint relaxes on the range of f2 m allowed for mod- 
els with w closer to —1, including ACDM. For the sake 
of argument, the 'best fit' models, where the likelihood is 
maximized, are also shown. We see more clearly that the 
shape test drives the preferred models away from w — — 1, 
towards w ~ —1/2. 

The concordance approach offers no such 'best fit', as 
it contains no procedure for weighting or combining data. 
However, in Figure |^ we carry out the exercise of artifi- 
cially shrinking all the error bars, to find the last remaining 
models. This is equivalent to imagining that all measure- 
ments have accurately determined the intended quantity, 
but overstated their uncertainties. This procedure nar- 
rows down to the same set of models, tt m ~ 0.33, but is 
not driven as strongly by the shape test as is the maximum 
likelihood procedure. 

The most potent of the intermediate red shift con- 
straints is due to type la supernovae, which we present 
in Figure [h]. In addition to the SCP results, the HZS 
group has presented two different analyses of their catalog 
of SNe, based on multi-color light curve shapes (MLCS) 
and template fitting; hence we show three SNe results. 
Carrying out a maximum likelihood analysis, all three give 
approximately the same result for the location of the 2a 
bound, favoring concordant models with low f2 m , and very 
negative w. It is interesting to observe that the SNe bound 
is consistent with the core of the low red shift concordance 
region, displayed earlier in Fi gure [^. Based on the SCP 
maximum likelihood analysis, Pcrlmutter et al. 1999 have 
reported a limit w < —0.6 at the la level. A \ z analysis 
of the same data gives a somewhat different result: the 
Fit C SCP data and the HZS data sets give comparable, 
although weaker, results to the likelihood analysis. In the 
spirit of conservativism, we have used the weakest bound 
which we can reasonably justify. Hence, for the concor- 
dance analysis, we use the 2a contour resulting from a % 2 
test. 

The statistical rate of gravitational lensing provides a 
counter to the trend towards low matter density. In Fig- 
ure |ll| we present the results of our analysis of the HST- 
SSS lenses. Other groups have come to similar conclusions, 
based on this and other l ens surveys . Our results are in 
excellent ag reement with Torres 



Waga 1996, Waga fc 



We have taken the attitude in our work that current 



1996 



observational uncertainties arc dominated by systematic 
errors, so that a conservative method of combining obser- 
vational constraints is by concordance. We apply the 2a 
limits for each individual observation to pare down the 
viable parameter range. However, it is interesting to com- 
pare this to what a naive maximum likelihood estimate 
(treating the errors as gaussian) would give. In Figure |^ 
we show the 2a contour in the fl m — w plane, where the 
remaining parameters have been marginalized. This pa- 
rameter region is only slightly smaller than that result- 
ing fr om concordance. It is reassuring that this technique 



ried out by Maoz fc Rix 1993 , Kochanek 199£ , Kochanek 
" |Falco et al. 1998| for ACDM. In general, there are 



yields approximately the same result, although one should 



fewer lenses observed than expected based on the volume 
red shift relation for a low density, A-dominated universe. 
The disparity between theory and observation is reduced 
as the matter density increases, or as w increases. 

We have evaluated the x-ray cluster evolution constraint 
using the observed abundance of rich clusters at z ~ 1. 
This test constrains the amplitude of mass fluctuations 
and the rate of perturbation growth. While it has been 
argued t hat this reduces to a bound on erg for ACDM 
models ( Bahcall fc Fan 1998| ), it has been show n that 
the bound depends on w for QCDM scenarios by Wang 
& Steinhardt 1998. In Figure O we show the conse- 



be cognizant of some of the pitfalls of both methods as 
discussed in the Appendix. In using the maximum like- 
lihood technique, we lose some of our ability to identify 



quence of the cluster evolution constraint on the concor- 
dance models in the Cl m — w plane; at this stage, the early 
formation of structure implied by the observations argues 
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against concordant quintessence models with an equation- 
of-state w J> —0.3. When the measurements comprising 
this constraint improve, we can expect a much more strin- 
gent result. Considering the hypothetical situation that 
future observations successfully reduce the systematic un- 
certainty to the present la level, the constraint boundary 
would shift to the small region with —0.8 <J w <; —0.5 and 
0.25 < fl m < 0.3. 

We have evaluated the high red shift constraint due to 
the select ensemble of C MB anisotropy m easureme nts, us- 
ing theCOBE, Python ([Piatt et al. 19971) MS AM flChcnp 
et a l l 1997|), QMAP flDevlin et al. 1993 , JHerbig et at 
1998f Qhyeira-Costact al. 1998|) Saskatoon flNcttcrficld"ct 
al. 19971), CAT QScott et al.l*99q ), and RING5M^citchci 
al. 19' j^) results. Based on a \ z test in 5T%, the set of con- 
cordant models projected down to the Cl m — h and fl m — w 
planes is unchanged from the low red shift concordance re- 
gion at even the la level. This "null" result from the CMB 
should not be too surprising; the current observational 
data is capable only of discerning a rise and fall in power in 
the Ci spectrum across I ~ 100 — 300. The results are un- 
changed if we include additional current CMB results, or 



use a x test m m (<^l )■> as suggested in Bond et al. 1998 



Rather, we must wait for near- future experiments which 
have greater ^-coverage, e.g. BOOMERANG, MAT, and 
MAXIMA, which are expected to significantly reduce the 
uncertainties. 

Since the submission of this manuscript, the data 
from the MAT (Miller et al, 1999) and BOOMERANG 
(Mauskopf et al, 1999) experiments have been released. 
However, neither significantly changes our results. 

Thus far we have applied the low red shift constraints 
in sequence with one of the other intermediate or high 
red shift constraints. It is straight forward to see how the 
combined set of constraints restrict the quintessence pa- 
rameter space. Taking the low red shift constraint region, 
which is shaped primarily by the BF, BBN, H, and as con- 
straints, the dominant bounds on the Q m — w plane are 
then due to SNe and lensing. The SNe drives the concor- 
dance region towards small f2 m and negative w; the lensing 
restricts low values of fi m . Putting these all together, an 
ultimate concordance test is presented in Figure 13. We 



see that the resulting concordance region in the Vt m — w 
plane is very similar to the core region obtained in Fig- 
ure |j| If the present observations are reliable, we may 
conclude that these models are the most viable among the 
class of cosmological scenarios considered herein. 

It is beyond the scope of the present work to deter- 
mine how well future observations will determine the val- 
ues of cosmological parameters in a QCDM scenario. How- 
ever, we are in a position to highlight those observations 
which appear well suited to testing the quintessence hy- 
pothesis. Clearly, the first g oal must be to distinguish 
Q from A (Huey et al. 1998). Observations which mea- 
sure the growth of structure at intermediate red shifts 
(z ~ 0.5 — 1.0) are best suited for the purpose. In this 
red shift regime, structure growth is still occurring for the 
A model, but has shut off significantly for the Q model. 
The Ly-a determination of the mass power spectrum am- 
plitude is at too high of a red shift to serve this purpose: at 
z ~ 2.5, evolution is still matter dominated for Q, m = 0.3 
models with w <. —1/3. Cluster abundances at z ~ 0.5 — 1 
and the supernovae magnitude-red shift relation are better 



suited to this goal. 

The shape of the mass power spectrum may prove to 
be a strong test of the QCDM scenario, if the obser- 
vation of a turnover in the power spectrum near k ~ 
0.02 — 0.0 6 h/Mpc and a break in the slope at higher wave 
numbers (Gaztanaga & Baugh 1998) bears out. It may 
prove difficult for the simplest, scalar field quintessence, 
or A for that matter, to generate su ch a feature. This is 
the subject of another investigation ( Zlatev et al. 1999 ). 

A test of the tracker quintessence scenario can be made 
by determining the change in the equation-of-state. If the 
equation-of-state can be measured at the present and at an 
earlier epoch, say z ~ 1, we can obtain a crude measure of 
the slope, dw/dt. Trackers have the special property that 
the equation-of-state becomes more negative at late times: 
w — > — 1 as D,q — > 1. A measurement of dw/dt > would 
argue against tracker quintessence. 

More exotic observational tests can be used to discover 
the presence of a quintessence field. For example, if the 
quintessence field is coupled to the pseudoscalar F I1V F IIV of 
electromagnetism a s suggested by some effective field the- 



ory considerations ( Carroll 1998 ), the polarization vector 



of a propagating photon will rotate by an angle Aa that 
is proportional to the change of the field value AQ along 
the path. CMB polarization maps can potentially mea sure 
the Aa from red shift ~ 1100 to now ( Luc et al. 1999 ) and 
distant radio galaxies and quasars can provide infor ma- 
tion of Aa from red shift a few to now ( Carroll 1998 ). If 
these two observations generate non-zero results, they can 
provide unique tests for quintessence and the tracker hy- 
pothesis, because tracker fields start rolling early (say, be- 
fore matter-radiation equality) whereas most non-tracking 
quintessence fields start rolling just recently (at red shift 
of a few) . 



5. CONCLUSIONS 

We have applied a battery of tests and constraints to 
the family of quintessence cosmological models, determin- 
ing the range of parameters which are concordant with 
observations. The most reliable constraints are those re- 
sulting from low red shift observations, and the COBE 
normalization of the mass power spectrum. These restrict 
QCDM models to a narrow range of parameters, char- 
acterized by low matter density, 0.2 < Q m < 0.5, and 
negative equation-of-state, — 1 < w <; —0.2. While the 
intermediate red shift results are still developing, the im- 
plications are very exciting. The SNe observations narrow 
the range of matter density near 17 m ~ 0.3 — 0.4, and force 
the equation-of-state to w <; —0.4. While this appears 
consistent with the core of the low red shift concordance, 
the potential for conflict is present if the matter power 
spectrum shape test demands w > — 1. Our results based 
on low red shift observations are given by Figures ||, ^[ 
adding the supernovae constraints, which are more recent 
and whose systematic errors have not been fully tested, 



produces the narrower range shown in Figure 13 



To what degree do current uncertainties in the Hub- 
ble parameter, the spectral tilt and other cosmic param- 
eters obstruct the resolution in wl To judge this issue, 
we have performed an exercise in which we fix h — 0.65, 
flf,h 2 = 0.019, and we choose the spectral tilt to insure 
that the central values of the COBE normalization and 
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the cluster abundance constraint are precisely satisfied. 
In Figures [li] and [lj, we show how different constraints 
restrict the parameter planes. Note first the long, white 
concordance region that remains in the £l m — w plane, 
which is only modestly shrunken compared to the concor- 
dance region obtained when current observational errors 
are included. The region encompasses both A and a sub- 
stantial range of quintessence. Hence, current uncertain- 
ties in other parameters are not critical to the uncertainty 
in w. The figure further shows how each individual con- 
straint acts to rule out regions of the plane. The color or 
numbers in each patch represent the number of constraints 
violated by models in that patch. It is clear that regions 
far from the concordance region are ruled out by many 
constraints. Both figures also show that the boundaries 
due to the constraints tend to run parallel to the bound- 
ary of the concordance region. Hence, shifts in the values 
or the uncertainties in these measurements are unlikely to 
resolve the uncertainty in w by ruling out one side or the 
other — either the constraints will remain as they are, in 
which case the entire concordance region is allowed, or the 
constraints will shift to rule out the entire region. 

New measurements not represented in this figure will 
be needed to distinguish A from quintessence. At this 
point, high precision measurements of the cosmic mi- 
crowave background anisotropy are the most promising. 

The results presented in this paper apply to quintessence 
models in which the equation-of-state is constant or slowly 
varying with time. In the latter case, setting w = w, given 
in equation ([?]), gives an excellent approximation to the 
observational predictions for a broad class of models. 

Particularly important classes of quintessence models 
are tracker and creeper fields. The tracker models are 
highly appealing theoretically because they avoid the 
ultra-fine tuning of initial conditions required by mod- 
els with a cosmological constant or other (non-tracking) 
quintessence models. An additional important feature of 
these models is that they predict a definite relationship be- 



tween the present day energy density and pressure, which 
yields a lower bound on t he constant, effective cq uation- 
of-state, near w ~ —0.8 ( Stcinhardt ct al. 1999| ). Note 
that the effective or averaged equation-of-state as com- 
puted from Eq. (m is about 10 per cent larger than the 
value of w today (given in Table I). In Figure [l6] we add 
this bound to the low red shift constraints, obtaining the 
concordance region for tracker quintessence. This region 
retains the core of our earlier low red shift concordance, 
and is consistent with the SNe constraints. Creeper fields 
occur for the same potentials as trackers, but the initial en- 
ergy density exceeds the radiation density at early times. 
The consequence is that the field rapidly rolls down the 
potential, towards a point which is only mildly (logarith- 
mically) sensitive to the initial conditions, where it sticks 
and is effectively frozen with constant potential energy at 
very early times. Hence, the creeper field has an equation- 
of-state w = — 1, and is effectively indistinguishable from 
a cosmological constant today. 

In Figure [l^ we combine all current observations on 
tracker models. Since these are arguably the best- 
motivated theoretically we identify from this restricted 
region a sampling of representative models, listed in Ta- 
ble I, with the most attractive region for quintessence mod- 
els being O m ks 0.33 ± 0.05, effective equation-of-state 
w k —0.65 ± 0.07 and h = 0.65 ± 0.10 and are consis- 
tent with spectral index n s = 1 indicated by the dark 
shaded region in Figure |l7|. These models represent the 
best targets for future analysis. The challenge is to prove 
or disprove the efficacy of these models and, if proven, to 
discriminate among them. 
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APPENDIX 

Throughout this paper, we have chosen to judge models by combining observational constraints according to the 
"concordance" method in addition to the maximum likelihood estimator (MLE) method. We have asserted that the 
concordance method is conservative, sometimes giving a more reliable judge of the situation than the MLE method, 
especially when the observational constraints may be dominated by systematic, nongaussian, and/or correlated errors. 
We advocate using both concordance and MLE methods, as we have done in this paper, and then analyzing the source of 
any discrepancy before determining which models should be ruled out. Since it has been commonplace to present MLE 
results alone, we thought it would be useful to illustrate some of the pitfalls that can arise. 

For this purpose, we employ a toy example in which we have two parameters, A and B, and two independent (obser- 
vational) constraints, represented as the 2a regions C\ and C2, which restrict the allowed ranges of parameters. This is 
meant to be a simplification of our real situation where we have a five-dimensional parameter space to analyze quintessence 
models, and we have many observational constraints. In our paper, we have tried to determine constraints in the Q m -w 
plane by projecting, effectively, from five-dimensions to two. In our toy model, we imagine projecting onto the A-axis 
to determine the constraint on A, indicated as a bar along the axis. Our interest is to compare the concordance region 
corresponding to 2a with the 95%CL contour from the MLE method. 

We present several simple examples in which there is a large disparity between the concordance and MLE procedures. 
In Figures |l9|, [l^, we represent a two dimensional parameter space, A-B, with two independent constraint regions, Cj.,2 
shown as shaded rectangles. The projection of the concordance and MLE regions, C conc and Cmle respectively, onto 
the horizontal axis are shown as thick strips. Although the following discussion is qualitative, the relative sizes of the 
projected strips are correct. 

Given the two observations, C\ and C2, the concordance region for A is obtained by: (a) finding their intersection in 
two-dimensions; and (b) projecting the two-dimensional intersection C conc onto the j4-axis to obtain a bar. Note that we 
do not project first, and then take the intersection. This method can lead to gross errors. For example, consider Figure 
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19, in which the C% and Ci have no intersection at all. By our method, the concordance region is properly identified as 
the null set, whereas projecting first and then finding the intersection would produce a considerable band of acceptance, 
a false conclusion. 

According to the MLE method, we are required to know the central value, //, of each region, and assume the errors, a, 
are Gaussian. In the following, we make the simplistic assumptions that the likelihood function for each observation is 
symmetric about the center of the constraint region, C\ or Ci. We weight each point x by a Gaussian /(x, jl, (?) for each 
of the two constraints, and identify Cmle, the contour of constant /1/2 or % 2 which contains 95% of the total probability, 
Jc fifadx = 0.95 J /1/2GS. This is a straightforward procedure. 

Case 1: Figure |l8| illustrates a case where C\ encompasses Ci- In this situation, the MLE method (indicated by 
the lower bar) produces a smaller acceptance region than the concordance method. If the errors are truly gaussian 
and uncorrelated, the conclusion based on MLE is the better representation of uncertainty. Note that the concordance 
region includes the MLE region plus an additional range of A; so, the error made with concordance is to include too 
much. However, no model is ruled out by concordance which intuition suggests ought not be eliminated. We call this 
"conservative," and our aim is to find a robust, conservative method. 

The MLE is not a robust, conservative method, as can be illustrated by the same figure. Suppose that the observations 
in C\ and C2 are suspected to be nongaussian or systematic or correlated. Then, we are clearly mistaken to rely on the 
MLE method and eliminate the range of A which lies within the concordance region but outside the MLE region. In 
this example, the difference is only modest, but if we were combining a number of observations, the MLE allowed region 
would be much tinier than the concordance region, and the error in trusting the MLE method would be more serious. 
This point is directly relevant to this paper. 

Case 2: The observational constraints C\ and C2 intersect in a small region or, as illustrated in Figure [l9[ have no 
intersection at all. This case is opposite to Case 1 in that the MLE method produces a larger acceptance region than the 
concordance method. For example, the concordance region in Figure [l9|is the null set, whereas the MLE contour suggests 
a large acceptance region. This is not a case where the MLE method is being conservative; rather, it is a case where 
the MLE method is misleading. Intuition dictates that the observational constraints are in conflict, and the concordance 
method reflects this conclusion by producing a null concordance region. The MLE method, taken at face value, suggests 
a broad range of agreement. To be fair to the MLE methodology, one is not supposed to accept the 95% MLE contour 
at face value. The contour represents a probability compared to the maximum likelihood point, and one is supposed to 
check that this point is indeed a good fit. In practice, though, this step is often ignored or discounted. For example, if 
the maximum likelihood point has a high value of x 2 by the conventional x 2 -test, this is often (properly) considered a 
problem due to underestimating experimental errors. But, as indicated by this example, the same statistical result can 
be an indication that there is a true contradiction between models and data, and that completely new models need to 
considered. Hence, a contradiction between concordance and MLE methods is a warning to examine closely the cause. 

Case 3: Suppose constraint C2 is obtained by combining many measurements with small statistical uncertainty but 
unknown correlated error. Then, in a MLE analysis, constraint C2 receives undue statistical weight. 

If we drop the simplistic assumption that the likelihood function for each individual constraint is symmetric about the 
center of the constraint region, which is rare, other kinds of discrepancies between concordance and maximum likelihood 
can occur. For example, the MLE region may be shifted with respect to the concordance region so that each test allows 
models which the other does not. 

As a prominent ex ample, the curren t , highly provocative me asurements of type la supernovae by the High-Z Supernova 



J1U HI |J1 VllllJlVjllV V^V UM. J..!. . U J. J.-.. ^■■l.lli.V. Il l,' ■ J. J. J. ^ J. LS._ y [J i. KJ V V/V ( t Ul V J. J. J. m k_l HI V111V11UU VI MJ UUJVVjI 11VVUIV J UJUUVy JL Xi^li I -J k_> LI yj IIVJ V t.lj 

Searc h Tea m (HZS: Riess et al. 199S, |Garnavich ct al. 1998 ) and the Supernova Cosmology Project (SCP: [Pcrlmuttcr 



et al. 1998) exemplify all three cases above. For the HZS data base d on MLCS (multi-color light curve shape) analysis 
or for SCP using the Fit C selected data set Pcrlmuttcr ct al. 199§| , the concordance region is significantly larger than 



the MLE acceptance region, as demonstrated by comparing Figures 10 and |13|. So, relying on the MLE 95%CL region 
eliminates models which formally pass the absolute x 2 test at the 95%level. While statisticians may argue that the MLE 
likelihood estimate is more reliable assuming uniform prior over the parameter space, in comparing qualitatively different 
models (which do not have uniform prior), the readers should beware that MLE can potentially rule out an entire model 
even if the model agrees at better than the 95% confidence level (as judged by ~x 2 ). Especially at this early stage when 
better data will soon be available, we advocate the more cautious, concordance approach. This first example is like Case 
1. We have also noted that the scatter in the supernova red shift - magnitude data is so wide that no model which we 
have tested passes a x 2 test with the SCP full data set. Hence, this is an example, like Case 2, where the concordance 
region is null but the MLE acceptance region is large. Finally, in a MLE or Fisher matrix analysis which combines the 
supernovae measurements with other observations, the fact that there are many individual supernovae with small reported 
uncertainty gives these measurements heavy statistical weight. However, as the survey teams admit, the measurement 
approach is new and there remains the possibility that as yet unidentified physical effects cause a systematic, apparent 
reddening of the data. As in Case 3, results obtained by simple statistical combination of supernovae data with other 
measurements should be viewed cautiously. 

Another feature of maximum likelihood analysis, well-known to practitioners but perhaps unappreciated by some, is 
the fact that the estimation of parameters in a multi-dimensional fit can be markedly different from the estimate when 
marginalizing over some param eters. For example, th e range of w spanned by the MLE 95% confidence region in the 
£l m -w plane (as considered by Pcrlmuttcr et al. 1999 and in this paper) is significantly narrower than the MLE 95% 
confidence region in the three-dimensional parameter space including h and significantly broader than the MLE 95% 
confidence region obtained by marginalizing over Q m and collapsing to a one-parameter fit to w. (This is due to the fact 
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that A\ 2 criterion for 95% confidence depends on the dimensionality of the parameter space.) 

These simple cases demonstrate the differences in the concordance and MLE procedures and, especially, some problems 
which can arise in MLE analysis. Because we maintain the position that systematic uncertainties dominate the errors in 
constructing the constraint regions, we advocate the concordance approach as being more conservative for cosmological 
analysis at the present time. In general, caution must be exercised when the two methods disagree significantly, and the 
source of the discrepancy must be understood in order to determine the true acceptance region. 
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Fig. 1. — The solid circles with lcr error bars are the Peacock APM data we use to test the power spectra shape. The five solid lines are 
quintessence models that pass the shape test with a confidence level of 95%. The model parameters for the black, blue, red, green, and purple 
curves are: w = -1, -1/2, -1/3, -1/6, 0; Qq = 0.70, 0.60, 0.55, 0.43, 0.20; T = 0.20, 0.26, 0.29, 0.37, 0.52. The bias factor is optimized for each 
model shown in order to obtain the minimum x '■ 

b = 1.01, 1.24, 1.46, 1.81, 2.49, respectively. All models have n s = 1 and Q b h 2 = 0.02. For 
comparison, the dashed line is a standard CDM model with best-fit bias 0.8 which fails the shape test. 
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Fig. 2. — The bulk velocity predictions as a function of radius (assuming a top hat window function) are shown for a representative set of 
QCDM models, and ACDM given at the end of this paper in Table I. Surrounding the best-fit quintessence model (Model 2), we have shown 
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Fig. 3. — The magnitude - red shift relationship determined by type la SNe is shown, for the HZS (using the MLCS analysis method) 
above, and SCP (full data set) below. The horizontal A(m — M) = reference line shows the prediction of an empty universe {Sltotal = 0), 
which has been subtracted from all data and theoretical curves. The thin dashed and dot-dashed curves show the predicted magnitude-red 
shift relationship for flat models with = 1 and S7 m = 1, respectively. The vertical offset of the data has been determined by minimizing 
the x 2 to the best fit ACDM model with Q, m = 0.3, which is given by the thick, solid curve. The predictions for quintessence (QCDM) models 
with w = —2/3, —1/3 for the same matter density are shown by the thick dashed and dot-dashed curves. 
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Fig. 4. — The upper panel compares the linear mass power spectrum at z = for the representative ACDM and QCDM models in Table I. 
All models are COBE normalized and satisfy the cluster abundance constraint on <rg. The solid and dashed curves have n s = 1; the dotted 
and dot-dashed curves have n B = 1.2. The shaded region in the top panel indicates where non-linear contributions are non-negligible. The 
lower panel shows the same power spectra projected back in time to red shift z = 2.5 and rescaled by the appropriate value of h at red shift 
z. Note that, in converting to z = 2.5, the abscissa in the lower panel has been rescaled so that it is expressed in terms of velocity; once this 
model-dependent rescaling is made, the models can be compared directly to the data. We show the constraints on the power spectrum, with 
lcr error bars, as deduced from the Ly-« forest. Among our representative models, the n s = 1 models are preferred over the n s = 1.2 models. 
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Fig. 5. — The bandpower averages used as the basis of the small angle CMB constraint are shown. For comparison, the black, blue, red, 
green, and purple curves are models 1-5, given at the end of this paper in Table I. The differences are small, but distinguishable in near-future 
experiments. All have acoustic peaks lying significantly above those of the standard cold dark matter model (thin black). The higher I results 
favor our representative models with n s = 1 over those with n s = 1.2. 




Fig. 6. — The projection of the concordance region on the f! m — h plane, on the basis of the low red shift observational constraints only, 
is shown. The observations which dominate the location of the boundary are labeled. 
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Fig. 7. — The projection of the concordance region on the Q m — w plane, on the basis of the low red shift and COBE observational 
constraints only, is shown. The observations which dominate the location of the boundary are labelled. If a wider range for the baryon 
density is allowed, such as 0.006 < Vl^h 2 < 0.022, the shape test and <rg constraint determine the location of the low Q m boundary, and the 
concordance region extends slightly as shown by the light dashed line. 
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Fig. 8. — The 2<r maximum likelihood contours in the Q m — w plane with the low red shift and COBE observational constraints only are 
shown. The dotted and dashed curves show the likelihood contours with and without the mass power spectrum shape test. The set of models 
which maximize the likelihood in each case are shown by the solid circle and the thick line. The shape test pushes models away from w = — 1, 
towards w ~ —1/2. The solid line shows the 2a allowed region according to the concordance region for comparison. 
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Fig. 9. — We carry out the exercise of shrinking the error bars on all measurements to obtain the equivalent 'best fit' models in the 
concordance approach. The surviving models have f! m ~ 0.33 and —0.9 < w < —0.7. This is similar to the best fit models obtained from the 
maximum likelihood method, however the concordance models are not as strongly affected by the shape test. 



No. 1, 2000 



COSMIC CONCORDANCE AND QUINTESSENCE 



23 



0.8 - 



0.6 - 



ft 



m 



0.4 - 



0.2 - 



0.0 



Supernovae Constraint 
based on 
Maximum Likelihood Estimate 




1.0 



■0.8 



■0.6 



-0.4 



■0.2 



equation-of-state w 



■o.o 



Fig. 10. — The 2cr maximum likelihood constraints on the Q m — w plane, due to the SCP (solid), HZS MLCS (short dashed), and HZS 
template fitting methods (dot-dashed). The light, dashed line shows the low red shift concordance region. 
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Fig. 11. — The 2a gravitational lensing constraint on the low red shift concordance region is shown. 
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Fig. 12. — The effect of the x-ray cluster abundance evolution constraint on the projection of the concordance region to the Q. m — w plane 
is shown by the dark shaded region. The light shaded region is due to the low red shift constraints only. 



2G 



WANG ET AL. 



Vol. 530 



0.8 



0.6 



ft 



m 



0.4 



0.2 



0.0 



LOW RED SHIFT + CMB 
+ SNe + LENSING 



H + BBN + BF 



lens//7 ( 




H+ c a 

BBN + BF 



1.0 



■0.8 



■0.6 



-0.4 



■0.2 



0.0 



equation-of-state w 



Fig. 13. — The dark shaded region is the projection of the concordance region on the Q m — w plane with the low, intermediate, and high 
red shift observational constraints. The dashed curve shows the 2a boundary as evaluated using maximum likelihood, which is the same as 
Figure 10 (See the Appendix for a comparison of the tests and a discussion of the pitfalls of the maximum likelihood approach.) 
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Fig. 14. — The concordance region (white) resulting if we artificially set h = 0.65 and £7(,/i 2 = 0.019 precisely and fix the spectral tilt 
to precisely match the central values of COBE normalization and cluster abundance measurements. The curves represent the constraints 
imposed by individual measurements. The curves divide the plane into patches which have been numbered (and colored) according to the 
number of constraints violated by models in that patch. 
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Fig. 15. — The concordance region (white) resulting if we artificially set Q(,h 2 = 0.019 and fix the spectral tilt to precisely match the 
central values of COBE normalization and cluster abundance measurements. The curves represent the constraints imposed by individual 
measurements. The curves divide the plane into patches which have been numbered (and colored) according to the number of constraints 
violated by models in that patch. 
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Fig. 16. — The concordance region based on COBE and low red shift tests for tracker quintessence is shown. The thin black swath along 
w = — 1 shows the allowed region for creeper quintessence and A. The equation-of-state is time- varying; the abscissa is the effective (average) 
w. 
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Fig. 17. — The overall concordance region based low, intermediate, and high red shift tests for tracker quintessence is shown. The thin 
black swath along w = — 1 shows the allowed region for creeper quintessence and A. The equation-of-state is time- varying; the abscissa is the 
effective (average) w as defined in Eq. (H). The dark shaded region corresponds to the most preferred region (the 2a maximum likelihood 
region consistent with the tracker constraint), f2 m 0.33 ± 0.05, effective equation-of-state w fs: —0.65 ± 0.10 and h = 0.65 ± 0.10 and are 
consistent with spectral index n = 1. The numbers refer to the representative models that appear in Table I and that are referenced frequently 
in the text. Model 1 is the best fit ACDM model and Model 2 is the best fit QCDM model. 
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Fig. 18. — The constraint regions C\ and C'2 do not intersect. In the concordance method for determining bounds on A, we first find 
the intersection of the Ci and C'2 in the full, higher-dimensional parameter space, and then we project that intersection region to obtain the 
constraint on A. In this case, the concordance region is null. The MLE method always allows some finite region of 95%CL. 
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Fig. 19. — The constraint regions C\ and C'2 overlap. The projection of the concordance and MLE regions onto parameter A, along the 
horizontal axis are shown by shaded strips. 
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Table I: A set of representative quintessence models satisfying all concordance constraints. For models #2 — 5, w is the effective equation-of-state. The age to is in Gyr, while H ta gives the age 
in units of the Hubble time, which is equivalently the instantaneous power at which the scale factor grows with time. The luminosity distance in units of the Hubble length is given by HqcLl at 
red shift z — 1; this quantity is directly related to the distance modulus for supernovae. For comparison with the open, empty reference cosmology used in Figure 3, HocIl — 1-5. The magnitude 
difference is then A(m — M) — 5 log 1Q (iJod£,/1.5). The rate of change of the equation-of-state for the tracker quintessence models, dw/dz, is evaluated at z — 0. The growth factor is defined 
as g = 5/d sc dm and the growth rate is / = d\nS/d\n a. The normalization of the mass power spectrum is given by 5h, where P(k) — 6^ H^™ 3 k ns T 2 (k)/(4ir) and T is the baryonic and dark 
matter transfer function. The rms mass fluctuation excess, <rg, o"5o, is evaluated for top hat window functions with radii 8, 50 Mpc/h. For the five representative models, <rg agrees with the cluster 
abundance constraint given in Eq. (5). Also, the values of Tjg = ^s^^ 6 as defined by Chiu et al. 1998 agree with their observational bound on rj based on cluster abundances and peculiar velocities. 
The rms bulk velocity, in km/s, is evaluated for top hat radii of R — 50, 150 Mpc/h, with 2a Maxwellian error bars. We have not listed any quintessence models located in the region excised 
between Figures 13 and 17 since models in this region are very similar to models 1 &; 2. For comparison, the last two entries show the model parameters and properties of an open CDM model, 
and SCDM. The tilt of the open model is not based on a particular open inflationary model but has been chosen so that the model is both cluster and COBE (following Bunn 8z White 1997) 
normalized. 



